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O ■ 

(*"*) . We make a new proposal to describe the very low temperature susceptibility of the doped Haldane 

^SJ ' gap compound Y2BaNii_iZna;05. We propose a new mean field model relevant for this compound. 

The ground state of this mean field model is unconventional because antiferromagnetism coexists 

f-s • with random dimers. We present new susceptibility experiments at very low temperature. We obtain 

a Curie- Weiss susceptibility x{T) ~ C /{Q + T) as expected for antiferromagnetic correlations but 

we do not obtain a direct signature of antiferromagnetic long range order. We explain how to 

\^ • obtain the "impurity" susceptibility Ximp(r) by subtracting the Haldane gap contribution to the 

total susceptibility. In the temperature range [1 K, 300 K] the experimental data are well fitted by 

Txi-m-piT) — Cimp (1 + Tlnip/r)"'^. lu the temperature range [100 mK, 1 K] the experimental data 

D ' are well fitted by Tximp(T') = j4In (T/Tc), where Tc increases with x. This fit suggests the existence 

}_( , of a finite Neel temperature which is however too small to be probed directly in our experiments. We 

also obtain a maximum in the temperature dependence of the ac-susceptibility x'iT) which suggests 

the existence of antiferromagnetic correlations at very low temperature. 
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I. INTRODUCTION 



Ch • • • • 

5^ . Disordered low dimensional spin systems are currently the focus of both experimental and theoretical interest. A 

^ ' very powerful theoretical technique used to describe these systems consists in iterating a real space renormalization 

group (RG) in which high energy degrees of freedom are progressively frozen out. The cluster RG was initially 

proposed by Dasgupta an Ma in the early 80's M and applied soon after by Bhatt and Lee to a three dimensional 

J> , (3D) model intended to describe Si:P H. In a series of articles, Fisher applied the cluster RG to disordered ID 

^^ '■ antifcrromagncts, including the random Ising chain in a transverse magnetic field and the random antiferromagnetic 

0^ , spin-1/2 chain ra. In these models, the cluster RG is so powerful that it allows to obtain exact results regarding 

^^ ■ two aspects of the problem: (i) the correlation functions associated with the approach of a random singlet fixed 

point; and (ii) the approach of a quantum critical point. It is fair to say that there exists now a detailed theoretical 

understanding of all random spin models in one dimension [0-w|, and even in higher dimensions [p[. Some models 

do not have a quantum critical point J 2 M while other models have a zero temperature phase transition which is 

controlled by the strength of disorder |^-3[. In particular, weak disorder is irrelevant in the random spin-1 chain but 

there is a quantum critical point at a critical disorder dc- For d > dc, the random spin-1 chain behaves like a random 

singlet ||5|,|6|, namely, the strength of disorder increases indefinitely as the temperature is scaled down, which is known 

— ^ . as an "infinite randomness" behavior [g|. 

r-| ' On the experimental side, it has become possible to fabricate low dimensional oxides in the last five years, and 

Q , dope them in a well controlled fashion. As a consequence, these systems give the unique opportunity of exploring the 

O ■ introduction of disorder in a spin gap state. It would be erroneous however to think that there exists a straightforward 

^ , relation between these experiments on quasi one dimensional oxides and the ID models that have been widely studied 

J^ ' by theoreticians in which the disorder is introduced in the form of random bonds. In fact, the ID models with random 

I ^S , bonds are not directly applicable to experiments. The reason is that the relevant models should incorporate the 

j^ ' following realistic constraints: 

(i) Doping in quasi one dimensional oxides is introduced by substitutions, and not in the form of random exchanges. 
The usual random bond models are not the most relevant ones. 



(ii) The temperature is finite in the experiments, even though often very small. 
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(iii) The realistic systems are not ID. Interchain couphngs become relevant at low temperature and can change the 
physics drastically. 

The problem is then to determine to what extent realistic models incorporating these three constraints still behave 
like the original ID random bond models. 

One of the recently discovered low dimensional oxides is CuGeOa ||] which has a spin-Peierls transition at Tsp — 
14 K below which the Cu02 chains dimerize with the appearance of a gap in the spin excitation spectrum. Soon after 
the discovery experimentalists began to study the effects of various substitutions on this inorganic compound. For 
instance, it is possible to substitute a small fraction of the Cu sites (being a spin- 1/2 ion) by Ni (being a spin-1 ion) JlG] 
or Co (being a spin-3/2 ion) ||Il|]. It is possible to substitute Cu by the non magnetic ions Zn ||l^-[l4[| or Mg [[l5[ . 
It is also possible to substitute some Ge sites (being outside the Cu02 chains) with Si fl^. The general feature 
emerging from the detailed experimental studies of the various substitutions is the existence of antiferromagnetism 
at low temperature. It has even been shown experimentally by Manabe et ai. that at low doping concentrations, the 
Neel temperature behaves like InT/v ^^ l/x |lj], therefore suggesting that there is no critical concentration associated 
with the onset of antiferromagnetism. 

The Haldane gap is another example of a spin gap in a low dimensional antiferromagnet [|l8| . Two inorganic 
Haldane gap antiferromagnets have been discovered recently: PbNi2V208 which has a spin gap A ~ 28 K [Q, and 
Y2BaNi05 which has a spin gap A ~ 100 K [|o). The substitution of the spin-1 Ni sites of PbNi2V208 with Mg (a 
non magnetic ion) generates long range antiferromagnetism at low temperature. The situation with Y2BaNi05 is not 
so well established. Previous works have failed to find long range antiferromagnetic order at low temperature [pl|~p4[. 
In particular /iSR experiments at 50 mK in Ref. |23] have reported a paramagnetic relaxation with a Mg doping of 
1.7% and 4.1%. We report here a study of the effect of Zn substitutions in Y2BaNi05. In agreement with Ref. p3[ , 
we do not find antiferromagnetic long range order down to 100 mK. But from the analysis of the experimental ac- 
and dc-susceptibility we deduce the existence of 3D antiferromagnetic correlations which rule out the "ID quantum 
criticality" scenario represented on Fig. 0-(a). 
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FIG. 1. Two possible scenarii for Y2BaNii_i:Zni:05. In (a) there are no interchain couplings {J± = 0) and the physics is 
governed by the approach to the quantum critical point of the disordered Haldane gap chain. A power-law susceptibility is 
expected in this situation at low temperature (Griffiths behavior) . In (b) interchain couplings play a relevant role and generate 
an unconventional phase in which antiferromagnetism coexists with random dimers. The approach to the Neel transition is 
controlled by the establishment of antiferromagnetic correlations. The purpose of the article is to show that the situation (b) 
can apply to Y2BaNii_a;Zni05. 



There have been several attempts to find a theoretical description of doped spin-Peierls systems. Fukuyama, 
Tanimoto and Saito have made the first proposal |E5[. Another route has been followed by Fabrizio, Melin and Souletie. 
In a series of article p^-|2^, these authors have provided a detailed scenario for disordered antiferromagnetism. 



Appealing aspects of their approach is that the doped spin-Peierls and the doped Haldane gap systems can be 
described in the same framework p8[ , and the model appears to be fully compatible with experiments. The relevant 
physics is determined by comparing two energy scales: 

(i) The coherence temperatm^e T* — Aexp [~l/(a::^)] which controls the formation of singlet correlations. 

(ii) The Stoner temperature Tstoncr ~ J±x^ which controls the formation of 3D antiferromagnetic correlations. 

In these expressions, x is the doping concentration, A is the spin gap, ^ is the correlation length of the pure gaped 
system, and J± is the interchain coupling. The specificity of Y2BaNii_j:Zn2,05 is that T* > Tgtoncr while T* <C 
Tstoncr in Cui_2:Zna;Ge03 and Pb(Nii_^Mgj,)V208 (see Rcfs. p^ , p8| |). As a consequence in Y2BaNii_2,Zn^05, 3D 
antiferromagnetic correlations coexist with non magnetic random dimer correlations (see Fig. n^(b)). Based on the 
calculation of the exchange coupling two spin-1/2 moments associated to the same impurity presented in Ref. [^, we 
propose a mean field model relevant for Y2BaNii_2:Zn2;05. From this mean field model we can obtain new insights in 
the nature of the low temperature phase, and compare the model to very low tempera ture susceptibility experiments. 
The article is organized as follows. The model is presented in section ||. Section |ll| is devoted to calculate the 
susceptibility and discuss the nature of the low temperature phase. In section ^ we present and analyze the very low 
temperature susceptibility experiments. Concluding remarks are given in section |y|. 

II. THE MODEL 
A. Defects in spin-Peierls and Haldane gap systems 
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FIG. 2. Schematic representation of two non magnetic impurities inserted in a spin-1 chain. The thick arrows are the spin-1 
moments. The thin arrows are the "edge" spin-1/2 moments. The exchange coupling the two spins Si and S2 is given by (|l|) 
and (y). We have represented an even segment. The exchange between Si and S2 is therefore antiferromagnetic. 

The model that should be used to describe doped Haldane gap compounds is well established. The introduction 
of a non-magnetic impurity in a spin-1 chain generates a pair of "edge" spin-1/2 moments on either side of the 
impurity |3C|-|3^ (see Fig. pi). The spin-1/2 moments associated with two neighboring impurities interact with the 
Hamiltonian 

H = J(/)Si.S2, (1) 

where the exchange is ferromagnetic for odd segments (Si and S2 are in the same sublattice), and antiferromagnetic 
for even segments (Si and S2 are in a different sublattice): 

J(O = (-)'Acxp(^-0, (2) 

where £, is the correlation length associated to the Haldane gap. We note J2 the ferromagnetic exchange coupling 
between two spins associated with the same impurity. Quantum chemistry calculations indicate that J2 is smaller 
than the interchain coupling J± in Y2BaNii_a;Zn^05 [p9J. The authors of Ref. p^] have obtained J± ~ 0.2 — 0.9 K 
and J2 < 0.1 K. This means that there is no temperature range in which the effective low energy model can be 
considered as one dimensional. On the contrary, when T is above J2 ^ J±_ the system can be represented by a model 
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FIG. 3. Schematic representation of two non magnetic impurities inserted in a spin-Peierls chain, (a) corresponds an even 
segment with two spin-1/2 moments. In this case the exchange (111), (Q) between Si and S2 is the same as for the Haldane gap 
chain, (b) corresponds to an odd segment with one spin-1/2 moment, (c) corresponds to an even segment with no spin-1/2 
moment . 

having J2 — Ja_ = 0. When T is below J2 ~ Ji_ the system develops directly three-dimensional correlations. As 
a consequence we will use a model with J2 = but with finite interchain correlations that will be treated in mean 
field. We note that the role of interchain correlations was already pointed out in Ref. |2J|. On the basis of DMRG 
calculations for open segments, the authors of Ref. |^^ could obtain a good agreement between the model and the 
experiments on Y2BaNii_a;Znj;05 above 4 K. The deviations occurring in the temperature range [2 K, 4 K] have been 
attributed to interchain correlations which play a central role in the model developed here. 

It will be fruitful to make a qualitative comparison with doped spin-Peierls systems, which are closely related to 
the doped Haldane gap systems p6[ . Non magnetic impurities introduced in a spin-Peierls system generate unpaired 
spin-1/2 moments (see Fig. S). Depending on the parity of the length of the segment and on the dimerization pattern, 
there can be zero, one or two spin-1/2 moments. If there are two spin-1/2 moments (see Fig. 0-(a)) the exchange (0) 
and (0) is the same as for the Haldane gap chain. To discuss the qualitative physics, we will make the approximation 
of including only the segments on Fig. 0-(a) that give rise to two spin-1/2 moments and neglect the other segments. 



B. Role of interchain interactions 

Interchain interactions can be described by introducing a correlation length S,i_ in the transverse direction. In two 
dimensions, the effective low energy Hamiltonian is a sum over all pairs of spins |p6[: 



where the exchange between the spin Si at coordinates (x.^, y^) and the spin Sj at coordinates (xj, yj) is given by 



(3) 




(4) 



where r^ — Xi — Xj and Vy — yi — yj. It has been shown from the RG approach in 2D that the Hamiltonian defined 
by (H) and (m in two dimensions has a Gnite randomness behavior [E6[. This indicates the existence of long range 
antiferromagnetism. In the present work, we do not address directly the model defined by (ph and (0) in 2D. Instead 



we propose that Y2BaNii_2:Zna;05 can be described by a mean field model in which we introduce a coupling to a 
molecular field —hg J2i ^i ■ The value of the molecular field is equal to the strength of antiferromagnetic correlations, 



which can be obtained from the susceptibility experiments in section IV. Since we use a mean field model, we cannot 
make the distinction between "antiferromagnetic long range order" and "antiferromagnetic correlations with a finite 
(but possibly large) correlation length" . For the sake of simplicity we assume that the staggered molecular field is 
uniform. But our approach can be extended to incorporate a distribution of molecular fields. The distribution of 
molecular fields can be calculated in a self-consistent way by imposing that the exchange field distribution is identical 
to the staggered moment distribution, up to a proportionality factor related to the strength of interchain correlations. 

C. Summary of the article 






FIG. 4. Schematic representation of the ground state obtained from the mean field approach. The even segments are coupled 
into non magnetic random dimers. The other spins give rise to antiferromagnetic long range order. 

Our strategy will be to solve the cluster RG in the presence of a molecular field for a model having J2 — (see 
Fig. y). Because this model reduces to a sum of independe nt two-s pin Hamiltonians, we call this model a "two-spin 
model". The Hamiltonian is given in Appendix ^ (see Eqs. Al, A5). 



The mean field model shows that antiferromagnetism coexists with random dimers at low temperature (see sec- 
tion |I|) which are due to the even segments coupled antiferromagnetically. This coexistence can be understood on 
a simple basis: at low temperature the even segments are coupled into non magnetic dimers while the odd segments 
are coupled into spin-1 objects. Interchain interactions generate antiferromagnetic correlations among the spin-1 
moments. The spin-1 moments can eventually order antiferromagnetically (see Fig. H). 

III. NATURE OF THE LOW TEMPERATURE PHASE: COEXISTENCE BETWEEN 
ANTIFERROMAGNETISM AND RANDOM DIMERS 

A. Susceptibility of the two-spin model 

The exact treatment of the two-spin model is briefly given in Appendix H. The total susceptibility is obtained by 
averaging Eq. A4 over all possible realizations of disorder, i.e. over all possible bond lengths: 

+00 
{{x{T)))=Y.V{l)x{l,T). 



1=1 



The bond length distribution is V{1) — x{l~xy^^ ~ xexp (— x/), with x the impurity concentration. The temperature 
dependence of the susceptibility is shown on Fig. g|. A maximum occurs in the susceptibility when the temperature 
is of the order of interchain interactions, below which the susceptibility decreases to zero. As it can be seen in Fig. |^, 
the contribution of even segments (being coupled antiferromagnetically) is much smaller than the contribution of odd 
segments (being coupled ferromagnetically). The reason is that even segments tend to form dimers which do not 
couple to a uniform magnetic field. Both types of segments couple to the staggered molecular field, which is why the 
susceptibility on Fig. ^ has a maximum when T ~ hs and decreases to zero at lower temperature. 

The temperature dependence of Tx{T) is shown on Fig. o for three doping concentrations corresponding to the 
experiments in section f^. At high temperature we recover the Curie constant of free spin-1/2 moments with a 
concentration 2x: 

At low temperature, Tx{T) is exponentially small because of the coupling to the staggered molecular field. 
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FIG. 5. Temperature dependence of the susceptibility of the two-spin model, with the parameters relevant for Y2BaNi05: 
^ = 6, A = 100 K, ha = 0.3 K, a; = 8%. We have also shown separately the contribution of even segments (coupled 
antiferromagnetically) and odd segments (coupled ferromagnetically) . 



0.04 
0.03 
^ 0.02 
0.01 h 



-I 1 I I I Mil 1 1 I I I III 

x=0.04 \ 

x=0.06 — X — 
h x=0.08 — ^ 




Temperature 

FIG. 6. Temperature dependence oiTxiT) for the two-spin model coupled to a molecular field. The parameters are identical 
to Fig. H. For each doping concentration we have shown by a dashed line the temperature dependence of the RG susceptibility 
( [L1[ ) which does not contain a coupling to the molecular field. 

B. Renormalization of the two-spin model 

Now we implement RG transformations for the two-spin model coupled to a molecular field h^ and having J2 ~ Q. 
At low temperature (T < hs) the spins are either frozen into dimers, or align antiferromagnetically. The coexis- 
tence between the two orders is therefore very natural in this framework. The detailed calculation in presented in 
Appendix ^. 

1, Order parameters 

The temperature dependence of the two order parameters is shown in Fig M for Y2BaNi05 and in Fig. for CuGeOa. 
It is clear that dimer formation is not important in CuGeOa while it plays a crucial role in Y2BaNi05. To check that 
this coexistence is not an artifact of the cluster RG, we have made a similar analysis with the direct diagonalization 
of the two-spin model (see Appendix H). Considering and eigenstate of the form a;|+, — ) + j/| — , +), we define 
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FIG. 7. Temperature dependence of the two order parameters (Paf{T) and ipdimiT) with the parameters relevant for 
Y2BaNi05 (see Fig. |5|). The energy scale corresponding to random dimer formation is larger than the e nergy scale at which 
antiferromagnetism appears. The calculation corresponds to the RG procedure presented in section [II B The parameters are 
given on Fig. H. 



<Pdim = -^{x-y) 



(6) 



as the order parameters associated to antiferromagnetism and dimer formation. In the presence of a pure dimer 
ordering, one has x = — y = l/\/2, leading to (/9af = and tpdim = 1- In the presence of a pure antiferromagnetic 
ordering, one has x = 1 and y — 0, leading to tpAF — 1 and (pdim = 1/2. The qualitative behavior is similar to the RG 
treatment (see Figs. P and Eol), therefore supporting our proposal that random dimers coexist with antiferromagnetism 
in the mean field model relevant for Y2BaNi05. 



2. RG Susceptibility 



The RG susceptibility in the temperature range 2hs < T < A is obtained as the sum of the contribution of the 
spin-1/2 and spin-1 moments that have not been decimated. We note ns(T) the density of spin-S* moments (number 
of spin-^* moments per unit length) . The total susceptibility reads 



X{T) = Y.^S{T) 



S{S + l) 
3T ' 



(7) 



As far as even segments are concerned, we find (see Appendix g) n^^^'^' (T) — x (T/A)^^, from what we obtain the 
susceptibility of even segments: 
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(8) 



As far as odd segments are concerned, we obtain from Appendix pt n~°i2 {T) = x{T/A)^^, and n]^ '{T) — 



1 - (T/A) 



x^ 



, from what we deduce 
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(10) 



The total susceptibility is Xtot{T) = XeveniT) + Xodd,i/2(r) + Xodd,i(^) 
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FIG. 8. Temperature dependence of the two order parameters ifiApiT) and ipdim{T) with the parameters relevant for Cu GeO ; 
(see Fig. |lO|). Dimer formation plays little role here. The calculation corresponds to the RG procedure presented in section [HE . 
The parameters are A — 44.7 K, ^ = 10, hs = 0.03 and x = 0.001. The values of hs and x are deduced from the low doping 
experiments in Ref. ]17|. 
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9. Antiferromagnetic and dimer order parameters with the parameters relevant for Y2BaNi05 (see Fig. m. The 



calculation corresponds to a direct diagonalization of the two-spin Hamiltonian, subject to a temperature-independent staggered 
molecular field. 
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(11) 



where the temperature T is such that 2hs < T < A. 

We have presented on Fig. pi a comparison between the direct diagonalizations and the RG susceptibility of the 
two-spin model. We see from this figure that the value of Tx{T) obtained from the cluster-RG is very close to the 
exact diagonalizations in the temperature range 2hs < T < A. We can make the following remarks: 

(i) In the absence of the staggered molecular field (i.e. in the absence of interchain couplings) Tx{T) tends to a 
constant when T ^ 0. This constant is equal to the Curie constant of spin-1 moments with a concentration 

x/2: 
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(ii) In the presence of a uniform staggered molecular field, the susceptibility is exponentially small at low tempera- 
ture: 
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FIG. 10. Antiferromagnetic and dimer order parameters with the parameters relevant for CuGeOa. The parameters are 
shown on Fig. H. 
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(13) 



if the applied magnetic field is parallel to the molecular field. If the applied magnetic field is perpendicular to 
the molecular field, the susceptibility tends to a constant: 



Tx{T) - XoT. 



(14) 



If the sample is a powder as in the experiments presented in section [V, the susceptibility given by Eqs.(|13| 
(M) should be averaged over all directions and we obtain 



rXav(r) r^ XxoT, 



(15) 



where A is smaller than unity. 



Therefore a relevant question for experiments is to determine whether Tx{T) tends to a constant (like in Eq. (|lz 
or tends to zero (like in Eq. (|l5|)) in the limit of zero temperature. From the answer to this question we can determine 
whether 3D antiferromagnetic correlations are important in Y2BaNii_j,Znj;05. 



IV. EXPERIMENTS 



We investigated the very low temperature susceptibility of the doped Haldane gap compound Y2BaNii_a;Zn2;05. 
We have performed ac- and dc-susceptibility measurements with a SQUID magnetometer equipped with a dilution 
refrigerator capable of measuring down to 100 mK and in magnetic fields up to 8 Tesla. The samples are the same as 
those used in Ref. |g^. In particular, there is one undoped sample with x ~ 0.00 and three samples with a Zn doping 
of x = 0.04, X = 0.06 and x = 0.08. 

Fig. 



11 



shows the temperature dependence of the dc-susceptibility in a log-log plot for the four samples where the 
low temperature data from 100 mK to 6 K are new, and the high temperature data from 2.4 K to 300 K correspond 
to Ref. [g^. Both measurements were taken in a static field of 0.1 Tesla. As can be seen, the data agree very well 
and no adjustments have been made in the overlap temperature interval. 

From Fig. y_2| showing l/x{T) versus T, we conclude that our experiments suggest the existence of antiferromagnetic 
correlations at very low temperature because the fits intercept the T axis at a negative temperature O. 

The susceptibility in the temperature interval [1 K, 10 K] can be more or less described by a Curie- Weiss law 
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From the insert of Fig. O showing Q versus x, we deduce the following: 
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TABLE I. Parameters of the best fits of the susceptibility data. O and C are the parameters of the Curie- Weiss fit ( |16[ ) in 
the temperature range [1 K, 10 K]. Timp and Cimp are the parameters of the power- law fit (|2d) of the "impurity" contribution 
to the susceptibility in the temperature ranee [1 K, 300 K]. Ai, Ci, A2 and C2 are the parameters of the two-exponential fit 
of the Haldane gap susceptibility (see Eq. (|l8|)). Tc and A are the parameters of the fit given by Eq. (Ell) in the temperature 
range [100 mK, 1 K]. 
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FIG. 11. Temperature dependence of the dc-susceptibility of the four samples {x = 0.00, x = 0.04, x — 0.06, x = 0.08) with 
a magnetic field H — 0.1 T. This is a log-log plot. 

(i) Q does not extrapolate to when x = 0.00. This is due to the fact that the pure sample contains already a 
certain level of disorder (see Ref. pJ). We can use x — 0.014 to describe the pure sample (see Fig. n3[). 

(ii) O increases with x which means that the strength of antiferroniagnetic correlations increases with the doping 
concentration. 

(iii) The 6 versus x variation is similar to Cui-a^ZuajGeOa (see Refs. |17| , p6| ). The similitude between Cui-xZnj^GeOa 
and Y2BaNii_a;Zna;05 is another indication in favor of the establishment of 3D antiferroniagnetic correlations 
in Y2BaNii_2:Zn2:05. 

Let us now improve the analysis of the low temperature susceptibility. We found that the experimental data in the 
temperature range [100 mK, 300 K] can be well described by a sum of two contributions (see Fig. nj): 



TX(T) = TXnnpiT) + rXHald(r), 



(17) 



where XHaid(2^) is independent of doping concentration and describes the susceptibility associated with the formation 
of the Haldane gap at high temperature. Ximp(T) describes the low temperature susceptibility associated with the 
introduction of doping. Following a suggestion by Souletie et ai. [p3 we have attempted to represent TxRaidiT) by 
the sum of two exponential contributions 



rxHaid(T) = Ciexp(-^) + C2exp(-^), 



(18) 



rather than the usual activated behavior 
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FIG. 12. 1/x versus T for x = 0.00, x = 0.04, x — 0.06 and x = 0.08. The dashed hne are hnear fits vahd in the temperature 
range [1 K, 10 K]. The fits intercept the T axis at a negative temperature which is a signature of antiferromagnetic correlations. 



S 

?5 




FIG. 13. Magnetic susceptibihty versus temperature at low temperature for x = 0.00, x = 0.04, x = 0.06 and x = 0.08. The 
curves have been fitted to the Curie- Weiss law ( |l6[ ) over the temperature range [1 K, 10 K]. The insert shows how O depends 
on the doping concentration x. To incorporate the existence of a finite level of disorder in the undoped sample, we made the 
replacement x = 0.00 —> x = 0.014. 



XHald 



(T)=Ar-i/2exp(_AHaid/T) 



(19) 



based on a quadratic expansion of the niagnon dispersion relation around the minimum at q ^ tt |35[ . The form (09^ of 
the susceptibility is a good description of Haldane gap chains at temperatures smaller than AHaid/2 [p6|J37[ . At higher 
temperatures, Souletie et al. ||3J] have found that (18) accounts extremely well for the results of exact diagonalizations 
and DMRG calculations on finite rings with an even number N of spin-1 Heisenberg spins, extrapolated to the limit 
N — > +00. The agreement between these numerical calculations and Y2BaNii„a;Zn2:05 is even quantitative since we 
found A2/A1 ~ 3.8 in Y2BaNii_:rZnx05, instead of A2/A1 ~ 3.86 from the numerical calculations, and C2/C1 ~ 0.25 
instead of C2/C1 ~ 0.32. Unlike Eq. dl9), the expression (O) tends to the Curie law in the high temperature limit. 
It is particularly interesting that the expression (nSh reaches a physically sound high temperature limit because this 
is in this limit that the Haldane term constitutes most of the magnetic signal. We took advantage of this situation 
to evaluate the Haldane gap term with the data above 1 K and subtracted the corresponding contribution to extract 
an improved estimate of the impurity contribution. The difference Ximp(r) — x{T) — XHaid(r) can be interpreted 
as the susceptibility of the spin-1/2 moments generated by the Zn ions. We have shown on Fig. WB the temperature 
dependence of Tximp{T). We verified that the high temperature limit of Tximp{T) corresponds approximately to a 
concentration 2x of spin-1/2 moments. 

To describe the variation of TximpiT) we use two different expressions in the two temperature intervals [100 mK, 
1 K] and [1 K, 300 K]. The impurity contribution shown in Fig. WH has been fitted by a power law in the temperature 
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FIG. 14. Fit of the experimental susceptibility in the temperature range [100 mK, 300 K]. We have shown Tx{T) versus T in 
a log-log plot. The susceptibility in the temperature range [100 mK, 1 K] has been fitted to an exponential term C3 exp ( — j?-) . 
The Haldane gap contribution ( jTsI ) is extremely weak at low temperature. For instance one has TxHaid(^) = 1-2 x 
10-6emu.K/Ba.mol if T = 10 K. 
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i-i 



10 

T(K) 

FIG. 15. Temperature dependence of rximp(r). The fits correspond to (EC 
concentration x. 




The insert show how 7 depends on the doping 



range [1 K, 300 K] : 



TXirapiT) = Cimp 1 



Ti, 



T 



(20) 



where Cimp and Timp are given in Table 



and 7 is shown on the insert in Fig. |15|. We have found for the temperature 
range [100 niK, 1 K] that the experimental data are better described by a logarithmic temperature dependence (see 
Fig. Ill): 



TXimp(T)=Alogio(^ 



(21) 



where A and T^ are given in Table |. Note that in this temperature range the contribution from the Haldane gap 
susceptibility (|l|) is negligible. 

From the two fits given b y ( |2(]| ) and ( |2l| ) we deduce that Tx{T) tends to zero in the limit T ^ 0. Regarding the 
discussion in section III B 2 we deduce that interchain interactions play an important role at low temperature. This 
validates the molecular field approach used in our theoretical description. Moreover the fit (po| ) suggests the existence 
of a finite Neel temperature Tc- We notice that 
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FIG. 16. Fit of rximp(r) to the form (|l|) in the temperature range [100 mK, 1 K]. 



(i) The value of Tc (see Table P is too small to be probed directly in experiments. 

(ii) Tc deduced from the fit ( pl| ) is one order of magnitude smaller than 8 deduced from the Curie- Weiss fit (O 
(iii) Tc increases with x, which is compatible with an interpretation in terms of a Neel temperature. 
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FIG. 17. Temperature dependence of the ac-susceptibility x'{T) with x = 0.08 for frequencies / = 0.11 Hz, / = 1.1 Hz, 
/ = 11.1 Hz. The temperature of the maximum of x'{T) is almost frequency-independent. 

Finally we also performed ac-susceptibility measurements which are shown on Fig. |l7|for three different frequencies. 
We have obtained a maximum in x'(T) at a temperature T ~ 160 mK. The temperature of the maximum in the real 
part of the susceptibility is frequency- independent. This indicates that in spite of a slowing down of the dynamics at 
very low temperature, Y2BaNii_a:Zna;05 is not a spin glass at very low temperature. This may be contrasted with the 
hole doped compound Y2-2;Ca2,BaNi05 JS^] (with < a; < 0.2) in which a spin glass thermodynamic transition has 
been reported at Tg = 2 — 3 K. The very low temperature ac-susceptibility of Y2BaNii_2;Znj;05 represented on Fig. |l^ 
is reminiscent of the dc-susceptibility of the mean field model (see Fig. 0) . We suggest that the very low temperature 
ac- and dc-susceptibilities (see Figs. M and O) can be interpreted in terms of antiferromagnetic correlations at low 
temperature. 
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V. CONCLUSION 

To conclude, we have presented a detailed study of antiferromagnetic correlations in Y2BaNii_2;Znj;05. Based on 
the quantum chemistry calculation in Ref. [^, we have proposed that Y2BaNii_j:Zn2:05 can be well described by a 
model in which interchain couplings are treated in mean field. The ground state of this model shows a coexistence 
between random dimcrs and antifcrromagnetism. In the presence of the staggered molecular field, the product Tx{T) 
goes to zero in the limit of zero temperature. 

We have presented very low susceptibility experiments and shown how to subtract from the total susceptibility the 
high energy contribution corresponding to the formation of the Haldane gap. The product rximp(T') tends to zero at 
low temperature. In view of the model, we conclude that this behavior is due to interchain correlations. At very low 
temperature the product Txii-ap{T) has a logarithmic temperature dependence in the temperature interval [100 mK, 
1 K]: Tximp(T) ~ Alog^Q [T/Tc). If we take for granted that this scaling is valid at lower temperatures we deduce 
the existence of a finite Neel transition temperature, which is in agreement with the assumption made in the mean 
field model. We have also obtained a maximum in the temperature dependence of the ac-susceptibility x'(T) which 
is an indication of the presence of antiferromagnetic correlations at very low temperature. 

An important question left open is to clarify the behavior of 2D and 3D models and determine whether there is an 
antiferromagnetic phase transition at a finite temperature. Using a mean field model we have shown here that there 
is a coexistence between dimer formation (being due mainly to the even segments coupled antiferromagnetically) and 
antifcrromagnetism (being due mainly to odd segments coupled ferromagnetically). If interchain interactions are very 
small, it is possible that singlets can be formed also among neighboring chains. These singlets would compete with 
3D antifcrromagnetism. An open question is to determine whether 3D antifcrromagnetism is stable against dimer 
formation even with very small interchain interactions. Answering this question may clarify why T^ (obtained from 
the fit ( pi| ) in the temperature range [100 mK, 1 K]) is one order of magnitude smaller than (obtained from the 
Curie- Weiss fit (|6|) in the temperature range [1 K, 10 K]). 

APPENDIX A: EXACT TREATMENT OF THE TWO-SPIN MODEL 

1. Even segments 

The Hamiltonian of even segments coupled antiferromagnetically reads 

n^J{l)^i.^2-hiSl + h2Sl (Al) 

where we included a staggered molecular field hs ~ {hi + ft.2)/2 and a uniform field h^ — [hi — /i2)/2. The uniform 
field is used to calculate the uniform susceptibility. There are two eigenstates [+,+) and | — ,— ) with an energy 
-E|_i_.+) = J/4 — hu, and i?|_,-) = J/4 + /i„. There are two other eigenstates |V'e) = Xe\+, — ) + 2/e| — , +) corresponding 
to e = ±1, with an energy 

J 



and 



E, = -j + -y.r- + 4h'„ 



l-,.^M^. m 



2 




Vi + 4W^' 



^ ^''^ Vl + 4/i2/j2 ^ ' 

The magnetization is given by the thermal average 

e-/3B|+,+> _e-/3^H,-> 
M(/,T) = ^^^ , (A4) 

where the partition function is Z{T) = Trexp (~f3Ti), with /? = 1/T the inverse temperature. 
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2. Odd segments 



The odd segments are coupled ferromagnetically: 

n = -I J(0|Si.S2 - hiSl - h2S^, 
hu in sublattice A and hi = /i2 = —hs + hu in sublattice B. The solution of (A5) is 



(A5) 



where hi ^ h2 — h 
straightforward. 



APPENDIX B: RENORMALIZATION OF THE TWO-SPIN MODEL 

We should distinguish between odd segments and even segments. We use the RG transformations given in Appen- 
dices |C| and 0. 

1. Even segments 



We start from the Hamiltonian given by Eq. Al in which the exchange is antiferromagnetic, and where we assume 
that there is no uniform magnetic field. There are three energy scales in the problem: 

(i) The temperature T. 

(ii) The exchange gap Aj = J = Aexp (— Z/C)- 

(iii) The staggered molecular field gap A^ — h^. 

There are two relevant length scales: 

(i) The thermal length ^t for which Aj = T. One has ^^ = Cln (A/T). 

(ii) The staggered molecular field length for which Aj — Ah. One has £,h = ^In {A/hg). 

We start the RG at high temperature, i.e. T ~ A, and decrease the temperature until ^t becomes equal to (,h- AH 
the pairs of spins such that I < ^t are frozen into dimers. Hence, the dimcr order parameter (pdim is given by 



V'dimCT) = £" ^V{l)dl = i 



-'^ 



x^- 



where T > h^ 



When the temperature becomes smaller than the staggered molecular field hg, the spins that have not yet been 
decimated (i.e. having I > ^h) contribute to the antiferromagnetic order parameter: 



<PAF(r) 



+ 00 



Vil)dl 



1 fh 



ih 



xi 



where T < hg 



2. Odd segments 



The odd segments are coupled ferromagnetically, with the Hamiltonian Eq. A5. At high temperature, the RG 



transforms the pairs of spin-1/2 moments having I < ^t into spin-1 moments, having a renormalized staggered 
molecular field equal to hg = {hi + /i2)/2 = /i^. The density of spin-1 moment reads 



fir I I 



-'? 



x^- 



where T > 2h^ 



These spin-1 moments have a gap Ah — 2hs. When T — 2hg, all of the spin-1 moments are frozen antiferromagnetically. 
The antiferromagnetic order parameter jumps from ipAwi'^hs -I- 0"*") = to 



ipAF{2hg-0+) 



2hs 

A 



< 
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When hs <T < 2hs, the pairs of spin-1/2 moments having a length ^2^^ < I < £,t are transformed into spin-1 objects, 
which are immediately frozen antiferromagnetically. The antiferromagnetic order parameter is 



VAf{T) = ifiAFi^K - 0+) + / " l:V{l)dl = i 



1 



x^- 



where h. < T < 21%. 



When T = hs, the survival spin-1/2 moments are directly frozen antiferromagnetically. The resulting antiferromag- 
netic order parameter is 



/■+°° 1 1 

ipAF(T) = ipAFihs) + / ^V{l)dl = -, where T < h^ 



APPENDIX C: RG TRANSFORMATIONS OF FERROMAGNETIC BONDS 



Let us consider the Hamiltonian 



n 



"•^1,281.82 — h\S^ — hiS2 



(CI) 



in which the spins 5*1 and 6*2 are coupled ferromagnetically. We replace the two spins S\ and ^2 by an effective spin 
5" = 5*1 -I- 5*2 with a Hamiltonian H. — —hS^, and need to calculate the renormalized field h. We use two different 
methods giving the same answer in the ferromagnetic case: (i) a semi-classical calculation using Holstein-Primakov 
bosons; (ii) a quantum mechanical calculation using Clebsch- Gordon coefficients. 



1. Holstein-Primakov bosons 



We replace the spin operators by their bosonic representation 5 f = 5*1 — a^a, S^ — y/2Sia, S^ — ^/2Sla^, and 
5| = 5*2 - b+b, S} = Vis^b, S2 = V^b+. The Hamiltonian Eq.^ becomes 



n = {-JS2 + hi) a+a + {-JSi + /i2) b+b + J^/S^ {ab+ + a+b) , (C2) 

where we discarded the classical energy contribution. The Hamiltonian Eq. |C2| is diagonalized by the unitary trans- 
formation a — c cos 6 + d sin 6, and b = —c sin 6 + d cos 6, where 



tan (26*) 



2J./S^ 



J {Si -S2) + hi~ h2 

Expanding the low energy Hamiltonian to lowest order in hi and /i2 leads to 

h{Si + S2) = hiSi + h2S2. 



(C3) 



2. Clebsch-Gordan coefficients 



We combine the two spins Si and ^2 to form a state of the type |5i, S'2|S'f , 5*1). The highest weight state is 
\S,S) = \Si,S2\Si,S2), and we need to determine the states 15, M) to calculate (5", M\ - hiSf - /i2 5'||S', M). 
It is straightforward to show that 



{S,S\ - hiSl - h2S'2\S,S) = -hiSl - h2S2. 



(C4) 



Identifying this with —hS, we obtain the renormalized magnetic field Eq. C3. One can check that Eq. C3 is also valid 
for spin-^* states having lower values of S^ . For instance, one finds 



\S, S-l) = J^\Si,S - 2\Si - 1, ^2) + J^\Si,S2\Si,S2 - 1) 



leading to 
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{S,S-l\- h^Sl - h2S'2\S, S-\) = -^^ (hiSi + h2S2) . 



Eq. C5 implies directly Eq. C3 
It can also be shown that 



\S,S - 2) - V Si2S ~ 1) I'^i''^^!'^! ~ -'-''^2) 



2W tttI^^^I^i, 52|5i - 1, ^2 - 1) 



5'(2S'-1)' 

/5'2(26'2-l) 
S'(2S'~1) 



5*1 , 5*2 1 5*1 , S'2 — 2) . 



We deduce 



(S*, S-2\- hiSl - h2S^\S, S-2) 



5-2 
S 



{hiSi + h2S2) , 



leading again to Eq. |C^. 

We are lead to conjecture that 



M 



{S, M\ - hiS'^ - h2S^\S, M) = -— {hiSi + h2S2) 



for any value of M. 



(C5) 



(C6) 

(C7) 
(C8) 

(C9) 



(CIO) 



APPENDIX D: RG TRANSFORMATIONS OF ANTIFERROMAGNETIC BONDS 



Let us consider the Hamiltonian 



Ti = ^1,281.82 — hiSl — /i2'S'2 



(Dl) 



in which the spins Si and S2 are coupled antiferromagnetically. We assume that Si > S2 and look for the effective 
Hamiltonian under the form —hS^, with S = Si — S2- 



1. Holstein-Primakov bosons 



To obtain the RG equation in the semiclassical limit, we represent the two spins by Holstein-Primokov bosons: 
SI ^ Si- a+a, St = ^/2S'la, S^ = ^/2S'la+ , and S^ = -S2 + b+b, S^ ^ V^b+, S^ = V^b. The Hamiltonian 
Eq. Dl becomes 



n = {JS2 + hi)a+a + {J Si - h2)b+b + JV5^(a6 + a+6+), 



(D2) 



where we retained only the quadratic contributions, and discarded the classical energy term. The bosonic Hamiltonian 
Eq. D2 is readily diagonalized by the Bogoliubov rotation a = ccosh6' + d"*" sinh6', with 



tanh(26l) = -- 



2J^/STS2 



J{Si+S2)+hi-h2 

Expanding the low energy Hamiltonian to lowest order in hi and /i2 leads to 

h{Si - S2) = hiSi - h2S2, 



(D3) 



which constitutes the classical renormalization equation. As we show in sections D2 and D2 , the quantum renormal- 
ization equation does not coincide exactly with the classical equation. 
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2. Clebsch-Gordan coefficients: coupling of a spin-s to a spin-1/2 

Considering 5i = s > 1/2 and ^2 ~ 1/2, we find 



|.-l/2,Af>^yi±|^^|.,l/2|M+ 1/2,-1/2) -y^^^^i^|.,l/2|Af- 1/2,1/2), 
from what we deduce (s - 1/2, M\ - hiSf - /i25'||s - 1/2, M) = -hM, with 



(D4) 



This RG equation is clearly not identical to Eq. D2 



3. Clebscli-Gordan coefficients: coupling of a spin-s to a spin-1 

Now let us consider Si — s > I and 5*2 = 1. Using standard orthonormalization and recursion techniques, we obtain 



\s-l,M) 



l{s- M){s- M + 1) 
2s(2s + l) 



s,l|M- 1,1) 



I {s + M){s - M) 
s(2s + l) 



s,l|Af,0) 



i{s + M+l){s + M) 



s,l|A/ + l,-l) 



y 2s(2s + 1) 

form what we deduce (s — 1, AI\ — hiSl — /i2>5'||s — 1, M) = —hM, with the rcnormalized magnetic field 

hs ~ hi{s + I) — ft.2- 



Comparing Eqs. D4 and Eq. D6, we are lead to conjecture that the general renormalization equation is 

HSi -S2 + l) = hiiSi + 1) - h2S2 (with Si > 52). 



(D5) 



(D6) 



(D7) 



The quantum RG equation Eq. |D7| docs not coincide exactly with the classical one (see Eq. DS). Nevertheless, the 
large-S* limit of the quantum RG equation Eq. D7 does coincide with the classical equation Eq. D3. 
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